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Abstract. 

With an increase in the application of ultra-wideband (UWB) electromagnetic 
pulses in the communications industry, radar, biotechnology and medicine, comes 
an interest in UWB exposure safety standards. Despite an increase of the scientific 
research on bioeffects of exposure to non-ionizing UWB pulses, characterization of 
those effects is far from complete. A numerical computational approach, such as a 
finite-difference time domain (FDTD) method, is required to visualize and understand 
the complexity of broadband electromagnetic interactions. The FDTD method has 
almost no limits in the description of the geometrical and dispersive properties of 
the simulated material, it is numerically robust and appropriate for current computer 
technology. In this paper, a complete calculation of exposure of the human eye to 
UWB electromagnetic pulses in the frequency range of 3.1-10.6, 22-29, and 57-64 GHz 
is performed. Computation in this frequency range required a geometrical resolution of 
the eye of 0.1mm and an arbitrary precision in the description of its dielectric properties 
in terms of the Debye model. New results show that the interaction of UWB pulses 
with the eye tissues exhibits the same properties as the interaction of the continuous 
electromagnetic waves (CW) with the frequencies from the pulse's frequency spectrum. 
It is also shown that under the same exposure conditions the exposure to UWB pulses 
is from one to many orders of magnitude safer than the exposure to CW. 
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1. Introduction 

While a large number of experiments have been performed in an attempt to provide an 
insight into the biological effects of the electromagnetic fields (Barnes and Greenebaum 
2006, Miller et al 2002, Polk and Postow 1995), the bioeffects of non-ionizing ultra- 
wideband (UWB) electromagnetic (EM) pulses have not been studied in as much 
detail as the effects of continuous-wave (CW) radiation. Nevertheless, as the interest 
for the application of UWB electromagnetic pulses increases, particularly in the 
communications industry and medicine, so does the interest to understand their 
bioeffects (Zastrow et al 2007, Ji et al 2006, Hu et al 2005, Schoenbach et al 2004). 
Only recently was the work in the UWB field in the United States declassified (Taylor 
1995), and only in 2002 were UWB pulses approved by the Federal Communications 
Commission in the U.S. for "applications such as radar imaging of objects buried under 
the ground or behind walls and short-range, high-speed data transmissions" (FCC 2002). 

A computational approach to the exposure of biological tissues to non-ionizing 
UWB radiation is more involved than the computation of exposure to CW radiation. 
In addition to a realistic description of the geometry, the physical properties of exposed 
biological material have to be known over a broad frequency range. The computation has 
to be numerically robust and appropriate for the computer technology of today. The 
finite difference-time domain (FDTD) method, introduced by Kane Yee in the 1960s 
(Yee 1966) and extensively developed in the 1990s (Sadiku 1992, Kunz and Luebbers 
1993, Sullivan 2000, Taflove and Hagness 2000), is a well known numerical method that 
satisfies these conditions. 

The response of a biological system to an EM pulse relates to the extent of 
conversion of EM pulse energy into mechanical or thermal energy as the pulse propagates 
through the material. Calculating the electric and magnetic fields around and inside 
the exposed material and combining these with the electromagnetic properties of the 
exposed material, the specific absorption rate, total deposited energy, and induced 
current inside the sample can be calculated. Finally, the interaction mechanism between 
the electromagnetic radiation and biological material can be modeled and understood. 
It is possible that the bioeffects of short EM pulses are qualitatively different from 
those of narrow-band radio frequencies. If, for example, the specific absorption rate 
standards defined for continuous radio frequency by the American National Standard 
Institute (ANSI 1992) or by the International Commission on Non-Ionizing Radiation 
Protection (ICNIRP 1998) are applied to the short UWB pulses, their meaning loses its 
clarity. Because of the pulses' short duration and as a result of the non-uniform power 
absorption, the deposition power and the induced current densities are large and can 
exceed locally even the highest allowed power limit of 4 W/kg (Simicevic 2007). At the 
same time, the energy from the exposure to an UWB pulse is too small to induce an 
increase in the temperature of the exposed sample due to, compared to CW, the much 
lower energy density of the UWB pulse. 

We have already applied the FDTD method to calculate the EM fields inside 
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biological samples exposed to nanopulses in a GTEM cell (Simicevic and Haynie 2005, 
Simicevic 2005). We included in the calculation the geometry and the physical properties 
of the samples to the fullest extent and applied a full three-dimensional computation, 
but have restricted the calculation to one sample (blood) and the frequency range up 



In this paper, the effects of the exposure of a human eye to UWB radiation in 
the frequency range of up to 90 GHz are calculated using the same three-dimensional 
FDTD computer code described and validated in the previous works (Simicevic and 
Haynie 2005, Simicevic 2005). The shape and the physical properties of the eye are 
described as accurately as allowed by the existing data. The computational space of 
30 mm x 26 mm x 26 mm was divided into ~ 2.0 x 10 7 cells and the computation was 
performed at the Louisiana Optical Initiative (LONI) cluster of supercomputers. The eye 
was exposed to the UWB pulses from three different frequency regions: 3.1 — 10.6 GHz 
(authorized for communications and radar imaging), 22 — 29 GHz (authorized for vehicle 
radars), and 57 — 64 GHz (authorized for unlicensed use) (FCC 2002). 

What follows is a detailed description of the computation: geometrical and physical 
properties of the eye, the properties of the UWB pulses, the details of the computation 
and data representation, and the results. 

2. The Shape and the Size of a Human Eye 

The shape and the size of a human eye, an example of which is shown in Figure [TJ 
varies from one individual to another and changes with age. While the theoretical eye 
models put the emphasis mostly on the optical properties of the eye (Norrby et al 2007, 
de Almeidal and Carvalho 2007, Siedlecki et al 2004, Lotmar 1971), they can as well be 
used in radiation protection (Charles and Brown 1975). 

For the application in the FDTD computing, the shape of the human eye was 




Figure 1. Schematics of the human eye. 



to 15 GHz. 
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Figure 2. Front of the human eye discretized by means of Yee cells. One cell in the 
figure represents 16 Yee cells used in the computation. 

discretized by means of Yee cells, cubes of an edge length Ax (Yee 1966). The two- 
dimensional shape of the eye shown in Figure [T] was normalized to the dimensions of 
the existing theoretical eye models, the pixels belonging to different eye tissues were 
retrieved by image processing software, and, because of rotational symmetry, rotation 
around the proper axis of rotation was performed to create the three dimensional model 
of the eye. To account for as small structures of the eye as, for example, the thickness of 
the cornea (only ~ 0.6—0.7 mm in size), and still be able to perform the computation in a 
reasonable time, the lengths of the Yee cell were chosen to be Ax = Ay = Az = 0.1 mm 
The computation was performed inside the computational space of 30 mm x 26 mm x 
26 mm in size discretized by ^ 2.0 x 10 7 cells. The front part of the eye constructed 
from Yee cells is shown in Figure [2} To make Figure [2] readable, 16 Yee cells used in the 
computation are presented as one cell in the figure. 

During the process of discretization of the human eye, the structures of the eye were 
separated into the tissues of known dielectric properties: sclera, vitreous humour, lens 
cortex, cornea, lens nucleus, muscles, and blood. The eye was embedded in the bone, 
muscle and skin tissue (wet and dry) . When dealing with UWB radiation it was crucial 
to have a proper description of the dielectric properties of the exposed tissues over the 
entire frequency range. They were taken from the established references (Gabriel 1996, 
Gabriel C et al 1996, Gabriel S et al 1996). 

From the relation between the size of a Yee cell Ax and the maximum frequency 
fmax allowed in the FDTD computation (Kunz and Luebbers 1993, Taflove and Hagness 
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2000) 



Ax^—^ , (1) 

±KJ J max 



where v is the speed of light in the material, the Yee cube edge length of Ax = 0.1 mm 
may be used to model exposure to the wave frequency of up to 90 GHz. This value is 
far greater than the frequency range of the UWB pulses used in the calculation. We 
have achieved an optimal agreement between geometrical description of the eye and the 
computational requirements of the FDTD method. 



3. Shape of the UWB Pulses 

The eye was exposed to vertically polarized EM pulses, each covering one of three 
different frequency regions: 3.1 — 10.6 GHz, 22 — 29 GHz, and 57 — 64 GHz. The shape 
of the pulse is obtained by the inverse Fourier transformation of the uniform spectral 
density of frequency width Af around the central frequency f c of the frequency region. 
In the time domain, this transform is described by the function 

„ sinh-A/Yt - t )l r xl (*-*q) 2 , x 

E = E 1 ° ;J cos[27r/ c (t - 1 , (2) 

nAf (t - 1 ) 

where E is the pulse amplitude, t is the time shift of the pulse, and a is the pulse 

(t-t ) 2 

width. The exponential term e ^ 2 allows for a smooth rise and fall of the pulse 
having, at the same time, a small effect on the pulse spectral density. The numerical 
value of the parameter t was the same for all the frequency regions, t = 3a. For the 
frequency region of 3.1 — 10.6 GHz: a = 166.7 ps, Af = 7.5 GHz and f c = 6.85 GHz, for 
the frequency region of 22 — 29 GHz: a = 66.7 ps, Af = 7 GHz and f c = 25.5 GHz; and 
for the frequency region of 57 — 64 GHz: a = 66.7 ps, Af = 7 GHz and f c = 60.5 GHz. 
The shapes and the frequency spectra of the pulses are plotted in Figure [3| 

Understanding the propagation of a UWB pulse is more demanding than 
understanding the propagation of a continuous EM wave. For some pulse frequencies 
the exposed sample looks like a conductor while for others it looks like a dielectric. The 
amount of the conduction currents compared to the displacement currents changes over 
the frequency spectrum of the EM pulse. In this paper we have restricted ourselves 
to studying the propagation of the pulses with almost uniform power spectra in the 
frequency range imposed by the FCC rules. 



4. Dielectric Properties of the Exposed Eye Tissue 

Electromagnetic properties of a tissue are normally expressed in terms of frequency- 
dependent dielectric properties and conductivity. When dealing with an UWB 
electromagnetic pulse, it is crucial that these properties are properly described over 
a large frequency range. One of the most accepted models describing the dielectric 
properties of a tissue is the Cole-Cole model. The references (Gabriel 1996, Gabriel C 
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Figure 3. Shapes and frequency spectra of the UWB pulses described by Equation [2j 
The amplitudes are normalized to 1 V/m. From the left are the shape and the 
spectrum for the pulse in the frequency region of 3.1 — 10.6 GHz, pulse in the region 
of 22 - 29 GHz, and pulse in the region of 54 - 67 GHz. 



et al 1996, Gabriel S et al 1996)) provide the four-term Cole-Cole parametrization of 
the dielectric properties of the ocular tissue used in the present work 

eiuo) = + V h ~r~ • (3) 

In this equation i = Soc is the permittivity in the terahertz frequency range, Asj, 

are the changes in permittivity in a specified frequency range, r k are the relaxation 
times, a is the ionic conductivity, and a k are the coefficients specific for the Cole-Cole 
model. They constitute up to 14 real parameters of a fitting procedure. 

Application of the Cole-Cole model is problematic for FDTD. It requires time 
consuming numerical integration techniques and makes computation unacceptably slow. 
If, instead of a Cole-Cole parametrization, the Debye model 

N A 

e(lj) =£oc + E 7—^ + " 4 

is used, the computation time can be reduced by an order of magnitude using very 
efficient piecewise-linear recursive convolution (PLRC) method (Luebbers et al 1990, 
1991, Luebbers and Hunsberger 1992, Kunz and Luebbers 1993, Taflove and Hagness 
2000). But, it is also important that the Debye model provides an equally accurate 
description of the physical properties of a biological tissue. 

In his thesis, G. R. Lugo (Lugo 2006) used an accurate, robust and efficient vector 
fitting technique (VECTFIT) (Gustavsen and Semlyen 1999) to replace the Cole-Cole 
parametrization by a multi-term Debye parametrization with no loss in the precision. 
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Table 1. Debye parameters for the tissues used in the computation. Static 
conductivity a is from Reference (Gabriel C et al 1996). 



His procedure is used in this paper. While the accuracy of the Debye parameterization 
increases with the number of terms, so does the computational time and the requirement 
on the memory, and an optimal number of terms has to be selected. It was found that 
in the frequency region considered in this paper the accuracy of the three-term Debye 
model 

£{UJ) =£oo + — — + — — + — — 5 

1 + lUJTi 1 + IUOT2 1 + MjOT 3 

is approximately the same as the accuracy of the corresponding Cole-Cole model. The 
parameters of the three-term Debye model for the tissue used in the calculation are 
presented in Table [T] and the comparison between the two models is shown in Figure [4} 



5. Field Calculation and Data Representation 

The quantities required to understanding the mechanisms and the consequences of the 
interaction between the electromagnetic radiation and biological tissue are obtained 
from the the values of electric and magnetic field components inside the tissue. Those 
components are accurate on the level of our knowledge of the physical properties of the 
exposed material. While the rapid developments in computer technology allow for the 
improvement of the geometrical description and for a longer exposure time calculation, 
for a reliable output, the knowledge of the physical properties of the material must be 
complete and correct. This is particularly important when simulating the interaction of 
UWB radiation with biological material. 

A qualitative understanding of the penetration of an EM pulse into the eye can 
already be achieved from the animated sequences of the pulse propagation. For the case 
of the propagation of the UWB pulses used in this paper, the complete animation can 
be accessed on-line (Simicevic 2007). Only selected snapshots of the intensity of the 
pulse penetrating the eye, shown in Figures [5| [6| and[7j can be presented in this paper. 
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Figure 4. Comparison of the relative permittivity of the tissues used in the 
computation parametrized by the Debye model (Equation |5j red solid line) and the 
Cole-Cole model (Equation [3j blue dashed line). Parametrization covers frequency 
range from 30 MHz to 100 GHz. 
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To understand the dosimetric implication from the exposure of a tissue to an 
electromagnetic pulse, or to model the interaction mechanism, a quantitative approach is 
needed. It mostly relates to the extent of conversion of EM pulse energy into mechanical 
or thermal energy inside the tissue. A detailed study is presented in the next section. 

6. Computation of Energy Dissipation for an Electromagnetic Pulse 

If an electromagnetic pulse propagates through biological material, the energy will be 
dissipated inside the material. The conversion of electromagnetic energy into mechanical 
or thermal energy inside a volume V is computed using equation (Jackson 1999) 

P= I J-EdV. (6) 

j v 

P is deposited energy per unit time, and J and E are, respectively, the current density 
and electric field inside the material. In the case of dispersive media, like the eye tissue, 
using Equation [6] is not trivial. First, the equation assumes a continuous distribution 
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Figure 6. Four different stages of penetration of EM pulse, in the frequency region of 
22 — 29 GHz, into the eye. The contours represent the instantaneous power density in 
decibels. Brighter regions correspond to higher power density. 



of charges and currents, which is not the case for the EM pulse propagation through 
tissue. Second, the current J consists of the conduction or ionic currents and of the 
displacement currents, which have to be correctly resolved (Jackson 1999). 

To properly calculate the amount of the absorbed energy inside the eye tissue, one 
first eliminates the current J in Equation ^ using the Ampere-Maxwell law (Jackson 
1999). The deposited energy per unit time is then expressed as a function of the EM 
fields inside the material 



dV. 



(7) 



D in this formula represents the electric displacement and B the magnetic induction. 
Volume integration covers the volume of the exposed object. 

Equation y consists of the term E x H , representing the energy flux density, and the 



term E • ^ 



H 



dB 
dt ' 



representing the total electromagnetic energy density. In the case 



of linear lossless non-dispersive media, the second term is interpreted as the difference 
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Figure 7. Four different stages of penetration of EM pulse, in the frequency region of 
57 — 64 GHz, into the eye. The contours represent the instantaneous power density in 
decibels. Brighter regions correspond to higher power density. 

of internal energy per unit volume with or without the presence of the EM field. In the 
case of linear dispersive media with losses, the case in this paper, such an interpretation 
is not trivial (Landau and Lifshitz 1960). In addition, the dispersion causes temporally 
non-local connection between D and E 

D(r, t) = s^Eir, t) + e f X (r)E(r, t - r) dr. (8) 

Jo 

Here e is permittivity of the free space, the permittivity at infinite frequency, and 
x( T ) is the susceptibility (Fourier transform of complex relative permittivity). This 
makes a direct calculation of total internal electromagnetic energy very difficult. 

In this paper, the dissipation of electromagnetic energy was calculated using the 
energy flux density. It was shown (Landau and Lifshitz 1960) that the formula for the 
energy flux density or Poynting vector S 

S = ExH, (9) 

is valid for variable fields even if dispersion is present. It is obvious that the difference 
of the energy flux entering a volume V and the flux exiting the same volume represents 
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the amount of the energy dissipated inside the volume. In the case of UWB radiation, 
this energy is obtained by integrating the Poynting vector over the pulse duration T 
and the impact area A enclosing the volume 



E 



It I a 



S • da 



dt. 



(10) 



Due to the absorption, this energy is ultimately converted into heat (Landau and Lifshitz 
1960). As a basic volume of interest we have chosen the Yee cell and performed numerical 



integration of Equation 10 



as 



E 



Pulse duration 



S-Aa 



Yee cell area 



At. 



(11) 



Aa is an element of the Yee cell area and At is the time step in the FDTD computation. 
Due to the offset in the location of the electric and magnetic field components relative to 
each other, and relative to the position of the Poynting vector, the value of the Poynting 
vector is obtained by using the average values of the adjacent field components. As an 
example, with the notation from Figure [8j the S x component of the Poynting vector S 
at a position (/, J, K) and time t is calculated as 

s.Ai,J,K) = E * {, ' J ' K) + *" i, - J - K+1> 

H,(I, J, K) + H Z (I, J, K+l) + HJJ-h J, K) + HJJ-h J, K+l) 

X 

4 

E Z (I,J,K) + E Z (I,J+1,K) 
2 

H y (I, J, K) + H y (I, J+l, K) + H y (I-l, J, K) + H y {l~\, J+l, K) 
x 4 W 

Components S y and S z are calculated in a similar way. 

The numerical accuracy of Equation [TT] was tested for the case of an EM pulse 
propagating through empty space where the deposited energy inside the Yee cell should 
be 0. In the worst case, the accuracy of the method, measured as the deviation from 
0, was found to be less than 0.02 %, a negligible error compared, for example, to the 
accuracy of the dielectric properties of tissues. 

Numerical computation of the absorbed energy was also validated by comparing the 



results obtained using Equation [TT] with those obtained by analytical solution. For a 
normal incidence of an EM wave on a infinite conducting surface of a dielectric constant 
e = £r — iei, the transmitted energy flux or transmittance can be calculated as 

(y/E R ~ l) 2 + £l 



T = l 



(V^r + I) 2 + ei 



(13) 



The dielectric constant for a muscle at a frequency of 7 GHz, the mid frequency of 
3.1 — 10.6 GHz range pulse, is e = 46.9 — il6.6. For this dielectric constant, Equation 



13 gives for the transmittance a value of 35%. Equation 11 gives for the transmittance 
of a 3.1 — 10.6 GHz range UWB pulse 33%. This is as good as expected agreement 
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Figure 8. The placement and labeling of the EM fields and Poynting vector 
components in the Yee cell, adapted from Reference (Kunz and Luebbers 1993). 



taking into consideration that the absorbed energy is computed only to the depth of 
20 mm and that the transmittance of the UWB pulse was compared to transmittance 
of a plane wave having the pulse's mid frequency. 



7. Dissipation of Electromagnetic Pulse Energy Inside the Eye 



The penetration of the EM pulses into a human eye is shown in Figures [9] and 10 , As 
expected, the penetration depth decreases as the pulse's frequency spectrum increases. 
In the case of the pulse in the lowest frequency range, as shown in Figure [9| significant 
energy penetrates about 10 mm into the eye where it is, as a result of the complex eye 
structures, relatively nonhomogeneously absorbed. Surprisingly significant contribution 
to energy absorption modulation comes from, in computation, typically neglected 
structures as, for example, the iris. Also, unexpectedly, most of the energy in the 
cornea is absorbed next to the eye lids. The 22 — 29 GHz range UWB pulse significantly 



penetrates into the eye only 2 — 3 mm (Figure 10). Most of the energy is absorbed by 
the cornea with the energy absorption hot spots again immediately above or below the 
eye lids. The 57 — 64 GHz range UWB pulse is almost entirely and uniformly absorbed 
by the cornea. While the energy absorption data exist for the entire eye volume, only 
selected data can be presented in this paper. 

Most of the experimental work on the effects of eye exposure to electromagnetic 
radiation in the frequency range covered in this paper was performed using continuous 
waves. We will compare the results of the UWB pulse propagation with the results of the 
CW, restricting ourselves to the frequencies above 10 GHz. Above 10 GHz experimental 
results exist for the exposure to 35 GHz waves (Rosenthal et al 1977, Chalfin et al 
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Figure 9. Contour plot of the density of the energy absorption across the vertical 
cross section of the eye in percentage of the total available energy for the 3.1 — 10.6 GHz 
range UWB pulse. Distribution of the absorbed energy along the mid-line is shown 
in the top figure. From left to right separate regions corespond to cornea, anterior 
chamber, lens cortex, lens nucleus, lens cortex (again), and vitreous humor. 



2002) and 60 GHz waves (Kues et al 1999, Kojima et al 2005). A short but detailed 
summary of the experiments, including the experimental apparatus, organic material, 
cellular environment, and possible electromagnetic interaction mechanisms can be found 
in Miller et al (2002). 

Searching for any abnormal effects resulting from the eye exposure to UWB 
radiation, the computed distribution of the absorbed energy in the case of UWB 
exposure was compared to the energy distribution in the case of CW. As shown in Figure 



11, no difference in the amount of absorbed energy was found when the 57 — 64 GHz 



frequency range UWB pulse was compared to 60 GHz CW, or when the 22 — 29 GHz 
frequency range UWB pulse was compared to 22.5 GHz CW, both corresponding to 
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Figure 10. Contour plots of the density of the energy absorption across the vertical 
cross section of the eye in percentage of the total available energy for the 22 — 29 GHz 
(left) and 57 — 64 GHz (right) frequency range UWB pulses. Distributions of the 
absorbed energy along the mid-line are shown on the top. From left to right separate 
regions corespond to cornea and anterior chamber. 



mid- frequency of the pulses spectra. When 35 GHz CW was compared to a 22 — 29 GHz 
frequency range pulse (bottom of Figure 11), as expected, the pulse penetrated deeper 
into the eye due to its lower frequency spectrum. No extraordinary effects outside the 
consequences from the Fourier decomposition of the UWB pulse were observed. Two 
factors played a role in the energy dissipation of the UWB electromagnetic radiation: the 
frequency spectrum of the pulse and the dielectric properties of the exposed biological 
material at those frequencies. 

The consequences of this study is summarized in the Table [2j As was shown, the 
energy in the material exposed to UWB electromagnetic radiation was absorbed in the 
same way as in the material exposed to a CW in the coresponding pulse frequency 
spectrum. The absorbed energy in the tissues was, in both cases, and in the same way, 
proportional to the total available EM energy. 

The total available energy in the case of CW can be easily calculated from the 
Poynting vector 

S = C ~yEl (14) 
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Figure 11. Top: comparison of the absorbed energy after exposure to 57 — 64 GHz 
frequency range UWB pulse (red line) and after exposure to 60 GHz CW (blue circles) . 
Bottom: comparison of the absorbed energy after exposure to 22 — 29 GHz frequency 
range UWB pulse (red line) and after exposure to 22.5 GHz (filled blue circles) and 
35 GHz CW (hollow blue circles). From left to right separate regions corespond to 
cornea, anterior chamber, lens cortex, and lens nucleus. 



where c is the speed of light, £q the vacuum permittivity, and E G the CW electric field 
amplitude. The energy of the UWB pulse can be easily obtained by integration 

S = ce J T E(t) 2 dt (15) 

over the pulse time duration T. E(t) is the electric field at the time t. It is easy to see 
that, for the same field amplitude, the energy carried by the short EM pulses depends 
on the pulse repetition rate and it is always much smaller than the energy carried by 
the CW. The results are summarized in Table [2| Assuming the same exposure time, 
for the pulse repetition rate of 1MHz, the energy absorbed by the exposed material 
will be r^j 0.01 % of the corresponding CW exposure, for the maximum repetition rate 
of 1 GHz, the energy absorbed will be ~ 10 % of the coresponding CW energy. Any 
possible absorbed-dose-related health effects as a results of the exposure to UWB will 
be therefore reduced by one to many orders of magnitude compared to the health effects 
caused by the exposure to CW. To quote from Miller et al (2002): " We want to close 
this section with a more general observation on the effects of UWB pulses. When we first 
started this work, there was speculation that even one UWB pulse might have serious 
effects on biological organisms. So far, we have exposed animals to over a quarter of a 
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Pulse Frequency Range (GHz) | RepRate: 1 MHz | RepRate: 1 GHz | 
3.1-10.6 0.011% 11% 

22 - 29 0.009 % 9 % 

57 - 64 0.009 % 9 % 

Table 2. Energy carried by the UWB pulse as a percentage of the energy carried by 
the CW of the same field amplitude, and in the same time interval. In one case 10 6 
and in the other 10 9 pulses were generated in one second. 



billion UWB pulses, and we have never seen any acute effects." 
8. Conclusion 

In this paper, we have performed full three dimensional FDTD calculations of the 
penetration of UWB electromagnetic pulses, authorized by the FCC for communications, 
radar imaging, and vehicle radars, into a human eye. Calculations included as detailed 
geometrical description of the eye as necessary and as accurate description of the physical 
properties of the eye tissue as allowed by existing data. The spatial resolution of 0.1 mm 
side length of the Yee cell allowed reliable calculation of up to ~ 90 GHz in frequency 
range in the dielectric. 

To minimize the computation time, the electromagnetic interaction with dielectric 
material was modeled using the Piecewise-Linear Recursive Convolution method 
(PLRC) for Debye media. Dielectric properties of the eye tissues in the frequency 
range < 100 GHz were formulated in terms of the Debye parametrization with the same 
accuracy as the accepted Cole-Cole parametrization. 

The energy absorbed after the UWB exposure was evaluated over the entire 
eye volume and compared to the energy absorbed after the exposure to continuous 
electromagnetic waves. The energy in the material exposed to UWB electromagnetic 
radiation was absorbed in the same way as in the material exposed to a CW in the 
coresponding pulse frequency spectrum. We have found that, assuming the same field 
amplitude and the same exposure time, any possible dose dependent health effects as 
a results of the exposure to UWB radiation will be reduced by one to many orders of 
magnitude compared to the health effects caused by the exposure to CW. This finding is 
in agreement with the experiments carried out at the US Air Force Research Laboratory 
at Brooks (Miller et al 2002). 

We can conclude that, based on the research described in this paper, any future 
applications of UWB pulses, for the purposes and in the frequency range approved by 
the FCC, pose less health risk compared to those applications being carried by the 
continuous electromagnetic radiation. 
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